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^ ; ABSTRACT 

fT^ ' High-precision astrometry can identify exoplanets and measure their orbits and masses, while coron- 

agraphic imaging enables detailed characterization of their physical properties and atmospheric com- 
positions through spectroscopy. In a previous paper, we showed that a diffractive pupil telescope 
(DPT) in space can enable sub-//as accuracy astrometric measurements from wide-field images by 

^^ ' creating faint but sharp diffraction spikes around the bright target star. The DPT allows simultane- 

5^ , ous astrometric measurement and coronagraphic imaging, and we discuss and quantify in this paper 

C^ ■ the scientific benefits of this combination for exoplanet science investigations: identification of exo- 

planets with increased sensitivity and robustness, and ability to measure planetary masses to high 
accuracy. We show how using both measurements to identify planets and measure their masses offers 
greater sensitivity and provides more reliable measurements than possible with separate missions, and 
therefore results in a large gain in mission efficiency. The combined measurements reliably identify 
potentially habitable planets in multiple systems with a few observations, while astrometry or imag- 
ing alone would require many measurements over a long time baseline. In addition, the combined 
measurement allows direct determination of stellar masses to percent-level accuracy, using planets as 
test particles. We also show that the DPT maintains the full sensitivity of the telescope for deep 
wide-field imaging, and is therefore compatible with simultaneous scientific observations unrelated 
to exoplanets. We conclude that astrometry, coronagraphy, and deep wide-field imaging can be per- 
formed simultaneously on a single telescope without significant negative impact on the performance 
of any of the three techniques. 

Subject headings: astrometry — Telescopes — techniques: high angular resolution — planets and 
satellites: detection 
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1. INTRODUCTION 

Among the existing techniques to identify exoplan- 
ets, astrometry and direct imaging are particularly well 
suited for identification and characterization of nearby 
habitable planets. Either technique can provide a full 
census of habitable planets around nearby (approxi- 
mately d < 10 pc) F,G,K main-sequence stars, provided 
that it achieves the required sensitivity (sub-/ias single 
measurement precision for astrometry, 10~^ raw contrast 
at w lOOmas for coronagraphy) . Extensive characteriza- 
tion of the exoplanet does however require both tech- 
niques, which provide complementary information. 

• Direct imaging with a high-contrast instrument 
is required to acquire spectra allowing character- 
ization of the planet's atmosphere. Direct imag- 
ing also reveals the exoplanet 's environment (other 
planets in the system, structure of the exozodiacal 
cloud), and allows measurement of the planet's ro- 
tation period and weather variability through time- 
series photometry. 

• AstrometrjQ is a promising technique to measure 
the mass of nearby exoplanets down to a fraction of 
an Earth mass, and this measurement is required 
to gain a physical understanding of the planet's 
surface and atmosphere. Although radial veloc- 
ity (RV) measurements can also measure planet 
masses, their ability to reach Earth-mass planets 
around Sun-like stars is limited by stellar jitter (we 
however note that an Earth-mass planet in the hab- 
itable zone of a low-mass star can produce a signal 
that exceeds stellar jitter). 

A planet's ability to retain an atmosphere and the 
composition of this atmosphere depend strongly on its 
mass. For low-mass planets, light observed by direct 
imaging may for example originate from the planet's sur- 
face, while planets with mass larger than a few Earth 
mass will likely retain a dense and opaque atmosphere 
hiding it. Meaningful interpretation of a planet's pho- 
tometry and color/spectrum will therefore require prior 
knowledge of its r nass, as colors alon e can be ambiguous 
(see, for example, iCahoy et al.[l201Cil) . 

Several astrometry ([Unwinl Hool iShao et al.l 120091: 



iMalbet et al.l 120 111) and direct imaging concepts 
jLcvinc et al 2009; Guvon et al. 2009; Trauac r et all 
[2010 : Shao fc Levine 2010: Clampin et al. 20061 have 
been proposed and studied for identification and charac- 
terization of nearby exoplanets, demonstrating the tech- 
nical feasibility of either measurement approach. While 
the scientific value of both astrometry and direct imaging 
for characterization of exoplanets is widely recognized, ei- 
ther path requires a dedicated mission with its own tech- 
nological challenges. Obtaining both astrometric and di- 
rect imaging measurements was until recently assumed 
to require two separate missions, and simultaneous de- 
velopment, construction and operation of both missions 
would likely be beyond available resources. 

^ In this paper, absolute astrometry (measurement of the abso- 
lute position of a star on the sky) is referred to as the astrometry 
measurement, while relative astrometry (measurement of the rela- 
tive offset between an exoplanet and its host star, usually from a 
coronagraphic image) is referred to as the coronagraphic measure- 
ment. 



I Guvon et al.l (|2012[ ) recently proposed using a diffrac- 
tive pupil telescope (DPT) to obtain both types of mea- 
surements simultaneously. The astrometric concept re- 
lies on simultaneous imaging of numerous background 
stars and diffraction spikes from the brighter target star 
on the same focal plane detector array. The diffraction 
spikes are created by a regular grid of dark spots on the 
telescope primary mirror. The submillimiter diameter 
dots cover a few percent of the primary mirror area, and 
their regular arrangement in the pupil plane creates in 
monochromatic light a regular grid of diffraction-limited 
spots in the focal plane. In broadband light, the spots are 
stretched into long thin spikes radiating from the central 
bright star (the target). Astrometric distortions in the 
imaging system (including the detector) are accurately 
measured by tracking the position of the spikes on the 
detector array, enabling high-precision astrometric cali- 
bration without requiring picometer-level stability of the 
telescope optics or focal plane detector array. The astro- 
metric measurement is performed by differential position 
measurements of the diffraction spikes originating from 
the bright central star targeted for planet search, relative 
to field stars. Because the spikes are generated by the 
primary mirror, such measurements are largely immune 
from astrometric distortions, as they equally affect the 
spikes and the field stars. The unique feature of DPT is 
that it allows simultaneous coronagraphic imaging within 
a narrow-field instrument along with astrometry and gen- 
eral astrophysics imaging with a wide-field camera. An 
optical design for the DPT is shown in Figure [T] for a 
1.4 m diam eter telescope . Thi s baseline configuration, 
adopted bv iGuvon et al.l (|2012D . offers a 0.2 ^as single- 
axis single-measurement astrometric precision for a 0.3 
deg^ field of view. This precision is a steep function of 
the telescope diameter, and improves, to a lesser extent, 
with the camera's field of view. 

While our first paper ([Guvon et al.ll201^ was focused 
on the technical aspects of the DPT concept and its ex- 
pected astrometric precision and accuracy, we discuss 
and quantify in this paper its ability to identify and 
characterize exoplanets and perform deep wide-field ob- 
servations that may not be related to exoplanet science. 
In Section [5J we quantify the required astrometric mea- 
surement precision for mass determination of exoplanets 
detectable by a coronagraph, and illustrate the scien- 
tific benefits of simultaneous coronagraphic imaging and 
astrometric measurement. We discuss in Section [3] the 
observation of planetary systems, which may consist of 
several planets. Finally we show in Section [H that the 
DPT concept allows deep wide-field imaging at the full 
sensitivity offered by the telescope, and that the impact 
of the added diffraction spikes on wide-field imaging sen- 
sitivity is negligible, even when the telescope is pointed 
at a bright star during exoplanet science observations, 
suggesting that the telescope could also be valuable for 
general astrophysics purposes. Detection and mass mea- 
surement of exoplanets with astrometry alone, as was 
previously considered for astrometry mission concepts, is 
briefly discussed in the Appendix and referred to when 
discussing the scientific merit of the proposed technique. 

2. EXOPLANET ORBITAL PARAMETERS AND MASS 
MEASUREMENT - SINGLE-PLANET CASE 

2.1. Science Goal and Representative Example 
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Fig. 1. — Conceptual optical design for the diffractive pupil telescope (DTP) proposed by IGuvon et ai.l II2012I ). The top part of this 
figure shows how light is shared between two instruments. The central field containing the bright star and its immediate surroundings is 
extracted at the telescope's intermediate focus and fed to a coronagraph instrument for high-contrast imaging. The wide well-corrected 
outer field is imaged onto a large focal plane detector array. Panels (a)-(d) show details of the wide-field image acquired in the final focal 
plane, (a) The wide-field image shows the diffraction spikes introduced by the dots on the primary mirror, (b) The central part of the 
target image, containing most of the flux, is missing from the wide-field image as it has been directed to the coronagraph instrument, (c) 
Faint diffraction spikes pave the rest of the field, (d) Faint background stars are imaged simultaneous ly with the diffracti on spikes. While 
images (a)-(c) arc simulated, image (d) was acquired in a laboratory demonstration of the technique IjBendek et aI.l[20111 V 



The science goals of the astromctric measurement ex- 
plored in this paper are twofold: 

• Assist the coronagraph to detect exoplanets. At 
the minimum, the astrometric measurement should 
confirm detections performed by the coronagraph 
and help constrain the orbital parameters of all 
planets identified by coronagraphic imaging. It 
may also uncover planets which are too close to the 
coronagraph's performance limits (especially inner 
working angle, IWA) to be firmly identified from 
coronagraph images alone, even if the planet ap- 
pears in at least one of the coronagraph observa- 
tions. 

• Measure the mass of all planets imaged by the coro- 
nagraph in the habitable zone of nearby stars 

In this section, we quantify the astrometric measurement 
precision required to meet these goals. We note that the 
second goal is more challenging than the first one, and 
that the performance required can therefore be derived 
from the second goal (mass measurement) alone. 

We assume that a 1.4 m telescope is used with a corona- 
graph offeri ng a 2 A/D IWA, follo wing the baseline design 
adopted bv IGuvon et al.l (|2012D . and inspired from the 
Pupil mappin g Exoplanet C oronagraph Observer mission 
concept (,Guvon et al.ll20ld[) . We consider that the num- 



ber and duration of observations is driven by the coro- 
nagraph instrument's goal to identify and acquire spec- 
tra of potentially habitable planets around a few high- 
priority targets: each observation is 48 hr long, and each 
high-priority target is observed approximately every two 
months. The expected astromet ric p recision of this sys- 
tem has been described in iGuyon et al. (2012), and was 
found to be 0.2 ^as per axis per measurement for a 0.3 
deg^ field of view camera. We note that the astrometric 
precision is independent of the coronagraph instrument 
design and performance, as it is achieved with an optical 
path which is separate from the coronagraph instrument. 
The relative mass estimation precision achieved for a 
fixed instrument design is a function of the planet type, 
stellar brightness (bright stars are easier, as their spikes 
are brighter in the astrometric camera), and star dis- 
tance (the astrometric signal is smaller for more distant 
stars). We choose to adopt a Sun analog at 6 pc from 
Earth as a representative example of a challenging tar- 
get for detection of an Earth-like planet with the coro- 
nagraph. The star apparent magnitude is my=3.7. We 
choose to place an Earth-mass planet at 1.2 AU around 
this star with a 1.3 rad inclination, therefore avoiding 
complications associated with the 1 yr period blind spot 
in the astromctric measurement. We use this example, 
for which the planetary system characteristics are listed 
in the top part of Table [TJ to quantify the relationship 
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between astrometric measurement precision and mass es- 
timation precision. 

2.2. Planet Mass Measurement 

In this section, we describe a simple model of simulta- 
neous observation of the system by the coronagraph and 
the astrometric camera. We then use the model to es- 
timate the astrometric precision required to achieve the 
science goals outlined in the previous section. 

Measured quantities. The planet orbit and mass are 
measured by combining the simulated astrometric and 
coronagraphic data. Coronagraphic and astrometric ob- 
servations are assumed to be simultaneous. For N 
epochs, the total set of measurements consists of AN 
values: 2N absolute astrometric measurements and 2N 
relative planet to star astrometric measurements derived 
from the coronagraphic images. Unless otherwise noted, 
simulations for this paper use a total number of 32 
epochs, regularly spaced every 57 days to span a total 
of 5 yr. With each observation assumed to be 2 day 
long, 28 targets could be observed at this cadence. 

Planet image position measurements from corona- 
graphic images. We assume that each coronagraph ob- 
servation yields the apparent position of the planet rel- 
ative to the star with a 2.5 mas Itr precision per axis, 
provided that the planet is outside the inner working of 
the coronagraph (assumed to be 130 mas). The assumed 
2.5 mas la precision per axis is equal to 1/30 X/D at 
A — 500/im, and would be achieved with an signal-to- 
noise ratio (S/N)=10 detection, or, in the photon-noise 
limit without background, with approximately 100 pho- 
tons. While in a 2 day exposure, the number of pho- 
tons collected from the planet is larger than 100, other 
sources of error will limit the measurement accuracy, 
such as centering error of the star on the coronagraph 
mask, and effect of uncorrected stellar speckles on the 
planet photocenter estimation. We have somewhat ar- 
bitrarily adopted 1/30 X/D as the single-axis measure- 
ment precision in this study, noting that a more de- 
tailed evaluation would be required to establish this mea- 
surement precision with high confidence. We note that 
centering errors in the coronagraph can be calibrated 
by introducing fai nt ghosts of the central s t ar, as pro- 
posed by [Si varam akrishnan fc Oppenheimeil ()2006[ ) and 
iMarois et al.i (|2006i ). and such a scheme is likely to be re- 
quired to reach the 1/30 X/D accuracy assumed here. If 
the planet is within the coronagraph mask, the position 
measurement error is increased to the size of the focal 
plane mask, essentially removing the contribution of this 
measurement from the overall solution. The apparent 
position of the planet on the sky is shown in Figure [2] for 
each of the 32 observations. The planet can only be seen 
by the coronagraph for 17 out of the 32 observations, as 
it is located within the coronagraph's IWA for 15 of the 
observations. 

Known variables. The following three parameters are 
assumed to be known: 

• Star apparent location on the sky (two variables). 
The parallax motion is therefore perfectly known 
if the star distance is known, and the aberration of 
light is therefore assumed to be perfectly known as 
it is only a function of source position on the sky 
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Fig. 2. — Coronagraphic (relative astrometric position of the 
planet referenced to the star) and astrometric (position of the star 
referenced to the sky) signatures are opposite in sign and of very 
different amplitude. The scale for the coronagraphic measurement 
is shown in the bottom and left, while the scale for the astromet- 
ric measurement is shown at the top and right. The dark shaded 
disk shows the zone within which the planet cannot be seen by the 
coronagraph. The dashed circle diameter is the one-dimensional 
standard deviation of a single astrometric measurement (0.2 ^as). 
The standard deviation of the coronagraph measurement (2.5 mas) 
is about half the size of the small green dots showing the planet 
position along its orbit. 



• Average radial velocity. The perspective accelera- 
tion effect produced by the relative radial motion 
of the star is therefore assumed to be known, and 
is not included in our simulations. For a star at 
6pc with a 20 km.s^^ velocity, the perspective ac- 
celeration ranges from to 0.2 mas.yr"^ depend- 
ing on the perspective angle (angle between three- 
dimensional velocity motion of the target and ob- 
server), and is maximum for a 45° perspective an- 
gle. If unknown, perspective acceleration will only 
affect the detection of planets with period longer 
than the observation time span, for which there 
is a degeneracy between the perspective accelera- 
tion and the exoplanet astrometric signature. For 
planets with orbital period comparable or less than 
the observation time span, the constant astromet- 
ric acceleration term is not strongly coupled with 
the planet signal, and can be included as a free pa- 
rameter to the astrometric data fit with little im- 
pact on planet detection and characterization per- 
formance. The degeneracy described above may be 
solved with rRV measurement of the host star to di- 
rectly measure the three-dimensional velocity and 
thus calibrate the perspective acceleration. With 
a 1 m.s— 1 RV measurement accuracy, and proper 
motion measurement obtained by astrometry with 
a 0.02 /las/yr accuracy as shown in Table [2l ana- 
lytical derivations show that the error on the per- 
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TABLE 1 
Observation Model 



Parameter 



Value 



Planetary system characteristics 

Star 

Distance 

Location 

Orbit semi-major axis 

Orbital period 

Planet mass 

Orbit inclination 

Orbit eccentricity 

Astrometric signal amplitude 



Orbit apparent semi-major axis 



Number of observations 

Coronagraphic image: planet position precision 
Coronagraphic image: inner working angle (IWA) 
Astrometry: single- measurement precision 



Sun analog 

6 pc 
Ecliptic pole 

1.2 AU 
1.3 yr 

1 Earth mass 

1.3 rad 
0.2 

0.5 lias 
200 mas 
Measurements 

32 (regularly spaced every 57 days) 

2.5 mas per axis (=3.6 mas in 2D) '^ 

130 mas 

Variable (driven by science requirement) 



"■ Corresponds to an S/N=10 detection with \/D = 80mas (single-axis astrometric precision 
is theoretically equal to {X/ D) / [it y' Npf^)) . For a photon-noise-limitcd measurement with no 
background, this would be achieved with Np^ =100 photons at 550nm for a 1.4m telescope. 



spective acceleration is 0.02 /las.yr^'^ (this error is 
maximum when the perspective angle is close to 
90°, corresponding to a near zero RV). Over a 5 
yr mission duration, the cumulative error around 
the mid-point observation reaches 0.06 /las at the 
beginning and end points of the observation series, 
and is therefore comparable to the astrometric ac- 
curacy assumed in this paper. An RV estimate 
with poorer accuracy is not very useful, as the per- 
spective acceleration is then better constrained by 
astrometry measurements alone. In this work, we 
do not take perspective acceleration into account 
in the global fit, implicitly assuming that (1) stel- 
lar RV is known to within 1 m.s— 1 or (2) the ob- 
servation time span is comparable or larger than 
the planet orbital period. Astrometric recovery of 
long-period planets (which is not the focus of this 
paper) would require a more careful treatment of 
perspective acceleration, and may be challenging 
as these planets are unlikely to be visible in coro- 
nagraphic images (reflected light decreases rapidly 
with orbital radius). 

In addition, it is assumed that observing epochs and 
standard deviation of all measurements (for both astrom- 
etry and coronagraphy) are known. 

Stellar mass. We assume that stellar mass is known 
to a 5% 1-cr accuracy prior to the observations, but we 
also treat it as a parameter to be solved for. This is done 
in our model by treating stellar mass as a free param- 
eter and adding, separately from astrometric and coro- 
nagraphic measurement, a single simulated stellar mass 
measurement with a 5% Ict accuracy. 

Correlation between measurements, systematic errors. 
We assume that all measurement errors are uncorrelated 
and all measurements are free of systematic errors. The 
validity of this assumption for the astrometric measure- 
ments is discu ssed in our previo us publication on the 
DPT concept (jGuvon et al.ll20l3 ) , and based on the fact 
that the astrometric measurement is inherently differen- 
tial. Existing performance estimates reflect our current 
assessment that dominant error sources in this system 



will be random, and that concurrent errors occurring in 
the different astrometric and coronagraphic fields of view 
will be primarily uncorrelated. However, we have begun 
to expand this assessment using detailed modeling and 
simulation, supported by lab testing, to characterize key 
error sources and quantify allocations of those errors in a 
comprehensive system error budget. We are particularly 
careful to identify any systemic error sources effecting 
the entire observational system, as well as the degree of 
correlation possible between errors both within our coro- 
nagraph and between our two concurrent observations. 

Free parameters to be solved for. A planetary model is 
constructed and linked to the observed quantities. The 
model is defined by 11 free parameters: central star dis- 
tance (one variable) , proper motion (two variables) . and 
mass (one variable), planet mass (one variable), and or- 
bital parameters (six variables). The values of the 11 free 
parameters are derived from the maximum likelihood so- 
lution given the noisy measurements listed above. 

Uncertainties for all 11 free parameters are determined 
using a Monte Carlo (MC) approach. For the single- 
planet case explored here, a total of 1000 simulations 
was used. Each data set is generated assuming a nor- 
mal distribution, centered on the perfect solution, with 
standard deviation equal to the respective uncertainties 
of coronagraphy and astrometry. 

2.3. Single- observation Astrometric Precision Required 
to Meet Science Goals 

To estimate the astrometric measurement precision re- 
quired to meet science requirements, simultaneous coro- 
nagraphic and astrometric measurements were simulated 
and fitted according to the scheme described in Section 
12.21 For each simulation, the single astrometric mea- 
surement precision was chosen in the range between 
(perfect astrometry) and 0.8 fias while the coronagraphic 
performance (both IWA and precision of the planet im- 
age position measurement) were kept constant. Figure [3] 
shows how the planet mass estimate changes as a func- 
tion of the level of astrometric error per measurement, 
and shows that a 10% relative precision in the estimate 
of the Earth-mass planet requires a 0.2 jias precision per 
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TABLE 2 

Uncertainties in the Combined (Astrometry and Coronagraphy) and 

Separate (Astrometry or Coronagraphy) Solutions'' 



Parameter 




1-a uncertainty 






Astrometry 


Astrometry 


Coronagraphy 




Only 


+Coronagraphy 


only"" 


parallax 


0.037 Atas 


0.035 Atas 


2.949 fias 


X proper motion 


0.017 ^as.yr-i 


0.012 /ias.yr-i 


1.304 /ias.yr-l 


y proper motion 


0.020 /ias.yr-i 


0.013 ^as.yr-i 


1.288 ^as.yr-l 


Planet mass 


0.132 MEarth 


0.098 Msarth 


5.355 MEarth 


Semi-major axis (SMA) 


0.0228 AU 


0.0052 AU 


0.0047 AU 


Orbital phase 


0.653 rad 


0.039 rad 


0.039 rad 


Orbit inclination 


0.0968 rad 


0.0065 rad 


0.0060 rad 


Position angle of SMA on sky 


0.111 rad 


0.0040 rad 


0.0039 rad 


Orbit ellipticity 


0.098 


0.0035 


0.0034 


Position angle of perihelion 


0.648 rad 


0.0034 rad 


0.0033 rad 


Stellar mass '^ 


0.05 Msun 


0.013 Msun 


0.012 Msun 



^ Assumed measurement precisions: 0.2 /las per axis per measurement for absolute 

astrometry, 2.5 mas per axis per measurement for relative astrometry derived from 

coronagraphic images. 

^ Assumes that astrometry is available at the 20fia.s per axis per observation from 

another mission to constrain parallax. Astrometric measurement at this precision level 

(2> fia.s) only affects estimates of parallax and proper motion, and has no significant 

effect on other measured parameters. 

•^ Assumed to be known to 5% accuracy independently of astrometry and coronagraphy 

measurements 
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Fig. 3. — Error on planet mass measurement as a function of 
single-measurement astrometric precision. The lower curve and 
points show the planet mass measurement precision when both 
coronagraphic and astrometric data are used to constrain the solu- 
tion. When only astrometric data are used, the mass measurement 
is less accurate. Each point in this figure shows the standard devia- 
tion of a set of 1000 solutions, each computed from an independent 
set of simulated data. 



astrometric measurement. 

The figure also shows that the mass estimate derived 
from the combined astrometry and coronagraphy mea- 
surements is more accurate than can be obtained from 
the astrometric measurements only. The difference be- 
tween astrometry alone and astrometry -I- coronagraphy 
is especially significant when high-precision mass mea- 
surement is to be obtained (left part of Figure [3]), as 
inferring the planet's mass from astrometry alone is then 
constrained by uncertainties in the stellar mass, assumed 
to be at the 5% level in this example. In this regime, 
improving the astrometric precision below 0.1 /xas does 
not significantly improve the planet mass measurement. 
With both astrometry and coronagraphy, the stellar mass 



is also measured, and high-precision planet mass mea- 
surement is possible. When the astrometric precision 
is relatively poor (right part of the figure) the astrome- 
try-|-coronagraphy measurements yield a planet mass es- 
timate standard deviation which is approximately half 
as large as would be obtained with astrometry alone. 
For example, with a 0.8 /ias precision per measurement, 
the planet mass is still constrained -^Q.'ihM Earth with 
the combined measurements, while astrometry only mea- 
surements would not be sufficient to unambiguously de- 
tect the planet and estimate its mass. The astrometry 
only points in the figure are consistent with the findings 
of the double blind study summarized in the Appendix: 
for the 0.8/ias precision per astrometric measurement, 
the mission astrometric S/N is ~ 4.5, insufficient (mis- 
sion S/N < 6) to unambiguously detect the Earth-like 
planet. This difference in mass measurement precision, 
due to the propagation of orbital parameters errors to 
the mass estimate, is described in the next section. In- 
terestingly, the relative difference in planet mass error 
between the two scenarios is the smallest when the sin- 
gle measurement astrometric precision is around 0.2 /xas, 
which is the case considered in this paper for numerical 
evaluations. 

Figure S] shows how the precision on the planet mass 
and semi major axis (SMA) varies as a function of both 
the relative and absolute astrometry measurement preci- 
sions. While the precision of the planet mass estimation 
(left) is mostly a function of absolute astrometry mea- 
surement precision, the precision of planet's SMA esti- 
mate is mainly a function of relative astrometric preci- 
sion. In both cases, both measurements however partici- 
pate to the parameter estimation. This figure illustrates 
the complementarity between the two measurements. 



2.4. How Coronagraphic Imaging Helps Astrometry to 
Measure Planet Mass 
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Fig. 4. — Error on planet mass measurement (left) and planet orbit semi major axis (right) as a function of single-measurement relative 
(vertical axis) and absolute (horizontal axis) astrometric precision. The maps contain white noise due to the finite number of simulated 
measurements used to produce the figure. 



While coronagraphic images (providing relative as- 
trometry between the star and the plan) alone do not 
constrain planet mass, they constrain the values for the 
orbital parameters much more precisely than possible 
with absolute astrometric measurements alone. The rel- 
ative astrometric signal given by the direct image has a 
large S/N advantage over absolute astrometry: for a 1 
MEarth planet on a 200 mas radius orbit, the assumed 
2.5 mas uncertainty corresponds to l/80th of the orbit 
radius. Constraining the orbital separation to the same 
level with astrometry alone would require a 0.006 /ias as- 
trometric precision (l/80th of the matching astrometric 
amplitude) . 

Solving for planet mass using the combined astrometry 
and coronagraphy measurements is therefore very pow- 
erful, and allows more accurate mass determination than 
would be possible with astrometry alone. 

• The coronagraphic images constrain the orbital pa- 
rameters and reduce error propagation from orbital 
parameter to mass estimate. For instance, [Torres! 
(119991 ) demonstrates how even a single direct image 
(providing an instant measurement of the projected 
orbital separation) dramatically improves the con- 
straint on the companion mass in binary systems 
for which only partial astrometric measurements 
are available. 

• With both coronagraphic imaging and astrometry, 
both the stellar and planet masses are directly mea- 
sured. Relative astrometry provides the total mass, 
while the absolute astrometric signal is driven by 
the mass ratio: both are required for accurate mea- 
surement of the planet's mass. 

The l-cr uncertainty for all parameters of the fit is 
shown in Table [2l for three cases: astrometric mea- 
surement only, combined astrometric and coronagraphic 



measurements, and coronagraphic measurement alone 
(assuming a lower precision astrometric measurement 
is available to constrain the parallax). Adding corona- 
graphic images to the astrometric data reduces the stan- 
dard deviation on orbital parameters by approximately 
a factor 10. The stellar mass is also directly measured, 
while with astrometry only, it is assumed to be known 
with a 5% standard deviation. The planet mass is esti- 
mated with a standard deviation below 0.1 MEarth with 
the combined data, while it would be 35% larger with 
astrometry only. Table [2] shows that the uncertainties 
in planet orbit geometry and stellar mass are identical 
(to the numerical accuracy of our test) in the astrom- 
etry -fcoronagraphy and coronagraphy only cases, illus- 
trating that these parameters are mostly constrained by 
the coronagraphic imaging data. 

2.5. Improving Mass Measurement Precision for 
Planets with a Nearly 1 yr Period 

The astrometric signature of a planet in a 1 yr period 
orbit can be partially absorbed in the parallax fitting of 
the astrometry measurements. With astrometric mea- 
surements only, the mass estimate error therefore grows 
as the planet period becomes closer to lyr for particu- 
lar orbit geometries. The problem is fundamentally that 
with finite S/N ratio the two signatures (parallax mo- 
tion and planet orbital motion) cannot be separated if 
their frequencies are close. The width of this blind spot 
is thus reduced as the astrometric measurements span 
a longer period, corresponding to a better resolution in 
the frequency domain. For a 5 yr mission duration, the 
astrometric blind spot covers approximately 5% of the 
habitable zone, so the detection rate is still good (approx- 
imately 95%) but is not 100%. The relative astrometry 
obtained from the coronagraphic images is not affected 
by parallax motion, and therefore yields an estimate of 
the orbital parameters that does not have any ambigu- 
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ity with parallax motion. With coronagraphic and as- 
trometric data, the blind spot is therefore greatly re- 
duced in size (although still not completely eliminated, 
as planet mass is still ambiguous if the planet orbit and 
phase matches exactly the parallax motion), especially 
for a short-duration mission. This benefit is critical for a 
combined astrometry + coronagraphy mission, as confir- 
mation that a planet is of sufhciently low mass to be po- 
tentially habitable needs to occur sufficiently early to al- 
low allocation of significant coronagraphic follow-up time 
(for spectroscopy and time monitoring). 

To illustrate how coronagraphic images mitigate the 
astrometric blind spot problem, we consider an Earth- 
mass planet at 1.01 AU from a Sun mass star (period = 
1.015 yr) at the ecliptic pole, orbiting a star at 6pc. We 
assume circular orbits for both the Earth and the tar- 
get planet, and a planet orbit phase equal to Earth orbit 
phase plus 1 rad, with a face-on orbit. This particularly 
unfavorable configuration is used as an example to illus- 
trate the blind spot mitigation thanks to coronagraphy, 
but is not a representative case of a habitable planet. The 
system is observed 32 times over 5 yr with observations 
regularly spaced in time. To demonstrate the robustness 
of the combined solution, we assume here measurement 
precisions that are poorer than in the previous section: 
standard deviation of 0.3 /ias and 5 mas per measure- 
ment per axis for respectively astrometry and coronagra- 
phy (instead of 0.2 /im and 2.5 mas, respectively). With 
only the astrometry measurements, the model described 
above yields a mass estimate with a standard deviation of 
4.17 Msarth'- the planet is not detected. With both the 
astrometric and coronagraphic imaging, the planet mass 
is estimated with a standard deviation of 0.16 MEarth- 

2.6. Discussion 

We have so far evaluated the scientific value of combin- 
ing astrometry and coronagraphy by considering a single 
example (5 yr observation of an Earth-like planet around 
a Sun- like star), and we now briefly discuss the applica- 
bility of the results presented in this section to other 
observation scenarios. While the orbit geometry (view- 
ing angle, orbit eccentricity) are not expected to signif- 
icantly impact detection sensitivity, the planet location 
(semi-major axis), planet and stellar masses, system dis- 
tance and observation time span must be considered to 
generalize the findings presented so far. 

The astrometric measurement is most sensitive to 
outer planets (the signal grows linearly with planet sep- 
aration and planet mass) provided that the observation 
time span is not significantly shorter than the orbital 
period. The coronagraphic observations are most sensi- 
tive to inner planets (reflected light signal proportional 
to the inverse square of planet separation) provided that 
the separation is larger than the coronagraph's IWA, and 
the measurement sensitivity is only weakly function of 
planet mass (the reflected light is proportional to the 
planet's squared radius, which goes as the l/3rd power 
of mass for low-mass planets, and increases less rapidly 
with mass for larger planets). The example considered 
in this section shows that combined astrometric -|- coro- 
nagraphic measurement can identify inner planets near 
the coronagraph's IWA and measure their mass to ap- 
proximately 0.1 Earth-mass accuracy. As the planet is 
moved outward, the astrometric sensitivity increases, al- 



lowing mass estimation to a « 0.12/SMA Earth mass, 
where SMA is the orbit semi major axis in AU. At a few 
AUs separation at most, the planet however becomes too 
faint to be imaged, regardless of its mass. The optimal 
mission duration is therefore close to the period of the 
outermost planets that can be imaged with the corona- 
graph, and there is little benefit in extending the mission 
duration further. 

Both astrometry and coronagraphy are more sensitive 
to potentially habitable planets as the stellar mass is re- 
duced (stronger astrometric signal, shorter period, and 
more favorable reflected light contrast). The reduced 
planet's angular separation however becomes challenging 
for the coronagraph, and the performance of the com- 
bined measurement scheme becomes driven by the coro- 
nagraph's IWA. 

As the system's distance to Earth is increased, both the 
astrometry and coronagraphic imaging signals become 
weaker, but this effect is much steeper for the corona- 
graphic imaging than for the astrometric measurement. 
The relative astrometric position derived from images 
suffers from the combined effects of smaller flux (planet 
flux scales as the inverse square of distance) and the an- 
gular resolution in absolute unit (scales linearly with dis- 
tance). The planet is also more likely to fall within the 
coronagraph's IWA, and the relative contribution of zo- 
diacal and exozodiacal backgrounds increases 

3. BENEFIT OF COMBINING ASTROMETRY AND DIRECT 
IMAGING FOR CHARACTERIZATION OF PLANETARY 
SYSTEMS 

3.1. Enhancement of Spectral Characterization 
Efficiency with Astrometry 

The benefits of astrometry for an exoplanet character- 
ization mission have been previously studied under the 
assumption that an astrometric measu rement would oc- 
cur prior to a direct imaging mission (jShao et al.l 120101 : 
lDavidsonll2011l : ISavranskv et al.ll2009H . and are twofold. 

• First, astrometric measurements would reveal the 
existence of planets and measure their masses, and 
therefore identify scientifically valuable targets for 
a future direct imaging mission aimed at acquiring 
spectra. 

• Second, the orbital parameters derived from the as- 
trometric observations could allow optimal timing 
of the direct imaging observations, when the planet 
is near maximum elongation. Telescope time can 
therefore be used more efficiently, as an imaging 
only mission can be rel atively ineffic ient to find 
planets, as quantified bv iBrownl (|2005H . 



iShao et aD ()2010t ) show that an astrometric precursor 
mission can reduce the coronagraphic observation time 
by a factor 2-5 thanks to these two advantages. The 
last advantage (ability to optimally time direct imag- 
ing observations) is however largely lost if several years 
separate the two missions, as the propagation of orbital 
parameters error over a long time span does not allow 
reliable prediction of the planet's apparent position. 

The first advantage (identification of scientifically in- 
teresting targets) is immune from error propagation of 
orbital parameters, and would allow considerable time 
saving, as the direct imaging mission would spend less 
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TABLE 3 

Multiple Planets System Characteristics 

Mass SMA Period Inclination 

{M Earth) (AU) (yr) (sin(j)) 



Planet 1 
Planet 2 
Planet 3 



1.0 
4.0 
16.0 



1.2 
1.8 
2.4 



1.31 
2.41 
3.72 



0.25 
0.25 
0.25 



time searching for planets and more time characterizing 
them. A direct imaging mission alone is not efBcicnt for 
identifying planets which spend a small fraction of their 
time outside the coronagraph's IWA, and the planet yield 
may be relatively small in that case. This is an unfortu- 
nate situation, as for most habitable planets the angular 
separation will statistically be close to the coronagraph's 
IWA. An imaging-only mission aimed at characterizing 
potentially habitable planets may also devote a signif- 
icant fraction of time acquiring spectra of planets out- 
side the mass range suitable for habitability, as apparent 
brightness is a poor proxy for mass. With prior astro- 
metric measurements, these risks are mitigated, and only 
a few direct imaging observations are required to recover 
the phase of the planet. Reconnecting the astrometric 
measurements with the direct imaging detection then al- 
lows accurate derivation of the orbital parameters if the 
time span between the measurements is large. 

In this paper, we propose to perform the astrometric 
and direct imaging measurement simultaneously, so we 
must reconsider the benefit of astrometry for improving 
the mission efficiency. Assuming a fixed mission dura- 
tion, identification of high-priority targets must be per- 
formed rapidly in order to allow sufficient time in the 
second part of the mission for extended characterization. 
An optimal observing strategy for a multi-year mission 
would be to first perform high-cadence, short-duration 
observations to identify potentially habitable planets, fol- 
lowed by several years of follow-up in order to character- 
ize the planet (s) through at least one full orbital period. 
Rapid identification of planets is especially challenging 
in multiple planets systems due to confusion between the 
astrometric signatures. In the following subsection, we 
use simulated observations of a multiple planets system 
to test how rapidly high-priority targets could be identi- 
fied. 

3.2. Multiple Planets 

We have simulated observations of a planetary system 
consisting of three planets around a Sun-like star at 6pc. 
The mass, SMA and orbital period for each of the planets 
are given in Table [3] 

Telescope and instrument design are identical to what 
was used in the previous simulations. Following the pre- 
cision requirement defined in Section [21 we assume a 0.2 
/ias uncertainty per astrometric measurement per axis. 

The goal of this section is to evaluate the combined 
measurement's scientific value when observing a realis- 
tic planetary system. We therefore include several plan- 
ets (three planets). While the astrometric precision is 
independent of the planetary system complexity, the 
relative astrometric measurement derived from corona- 
graphic images will be affected by confusion issues (such 
as the potential inability to link planets to their images 
for widely separated observations, relative astrometric 



error due to unknown planet affecting the planet image 
photocenter measurement, and spatial structure in the 
exozodiacal cloud). For simplicity, we have chosen to not 
include these effects in our model, but to arbitrarily ac- 
count for them by adopting a poorer relative astrometry 
precision. The relative astrometry measurement preci- 
sion is thus assumed to be 0.0088 arcsec (1/lOth of the 
diffraction limit at A — 550 nm) per axis per measure- 
ment for each of the three planets, provided that the 
planet image is outside the coronagraph's IWA. 

To illustrate the complementarity of astrometric and 
coronagraphic measurements, and the value of perform- 
ing both measurements, we have simulated two different 
mission scenarios. 

• Astrometric mission. One observation is performed 
every month, with an 0.2 /xas precision per axis. 

• Astrometry + coronagraphy mission. Every two 
months, an astrometric measurement at the 0.2 
/.ias precision per axis is performed, simultaneously 
with a coronagraphic observation. The planet po- 
sitions are measured to 0.0088 arcsec precision per 
axis only when outside the coronagraph's IWA. 

A number of simplifications to the observing model 
have been used for the simulations. 

• The orbits are assumed to be circular and coplanar. 
This leaves a total of 12 free parameters: mass, 
SMA and orbital phase for each of the three plan- 
ets, in addition to the stellar mass, its distance and 
the common orbital plane inclination. 

• The planets are assumed to recognizable in the 
coronagraphic images: provided that a planet is 
seen outside the coronagraph's IWA, it is assumed 
that it cannot be confused with other planets (this 
would be the case if each planet had recognizable 
colors) 

• The planets are detected in the coronagraphic im- 
ages if and only if outside the coronagraph's IWA, 
and the precision of the planet's measured image 
position is independent of its location in the im- 
age or phase. In reality, the measured planet po- 
sition would gradually become more uncertain as 
it is closer to the coronagraph's IWA, and con- 
trast/sensitivity limits would prevent its detection 
in coronagraphic images when illumination fraction 
is small. This assumption leads to an underesti- 
mate of the required number of observation epochs. 

We use again an MC approach to determine how each 
mission scenario would perform on the simulated three- 
planet system. A total of 200 MC simulations were con- 
ducted in each case, using a normal distribution to pro- 
duce each noisy data set. For each of the 200 MC simula- 
tions, a random initial solution is chosen with three plan- 
ets as a starting point for a x^ minima search, with ini- 
tial masses randomly chosen in the 0.5Mj-1.5Mi range, 
where Mi is the true mass of planet i, and with ini- 
tial SMA randomly chosen in the Q.bSMA, - 2.?)SMAi 
range, where SMAi is the true SMA for planet i. Simi- 
larly, the initial orbital phase is randomly chosen within 
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Fig. 5. — Result of the astrometry mission scenario, observing a three planet system around a Sun-like star at 6pc. Values to be 
measured are shown on the left with measurement error bars (Icr). Two hundred realizations of the observations are simulated, and the 
best three-planet solution is shown on the right for each simulation in the semi major axis— mass plane. Total mission duration is increased 
from top to bottom, from 1 yr to 4 yr. The three dark circles indicate the actual configuration for the three-planet system. The initial 
starting point is the same for all scenarios: the first six absolute astrometry values are identical in all cases. However, the astrometric 
values shown in the left panel are corrected for uniform drift (proper motion) during the measurement period, resulting in different shapes 
and amplitude for the part of the astrometry trajectory curves common to all scenarios (first year). 
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Fig. 6. — Result of the astrometry -|- coronagraphy mission scenario, observing a three planet system around a Sun-like star at 6pc. Values 
to be measured are shown on the left with measurement error bars (Icr). Two hundred realizations of the observations are simulated, and 
the best three-planet solution is shown on the right for each simulation in the semi major axis (x axis) - mass (y axis) plane. Total mission 
duration is increased from top to bottom, from 1 yr (top) to 4 yr (bottom). The gray disk indicates the blind area (within coronagraph 
inner working angle) for the planet-to-star relative astrometry measurements. The three dark circles indicate the actual configuration for 
the three-planet system. 
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1 rad of the true orbital phase. The sine of the plane- 
tary system inclination sin incl is randomly chosen with 
uniform probability distribution between and 1. This 
choice of initial conditions ensures that the search for 
the optimal solution is allowed to find minima outside 
the true solution of they exist. A simulated annealing is 
used to find the solution with the smallest x^j with 5e6 
iterations. The estimated values for the 12 free parame- 
ters are derived from this maximum likelihood solution. 
Results of this simulation are shown in Figures [S] and [5] 
for the astrometry-only and astrometry-fcoronagraphy 
mission scenarios. Measurements and solutions are 
shown for variable mission durations, and solutions are 
shown in the two-dimensional SMA versus planet mass 
space. 

3.3. Results 

Figure [5] shows how well the planetary system char- 
acteristics are estimated for the astrometric mission, as- 
suming a total mission duration ranging from 1 yr (top) 
to 4 yr (bottom). With a 1 yr long mission, none of 
the planets is unambiguously detected, and the measure- 
ments only show the existence of at least one planet with 
orbital period exceeding the mission duration. With a 2 
yr mission duration, the astrometric data can exclude the 
presence of massive planets at small angular separation: 
no solution places a planet of more than 10 Earth mass 
within the central 1.5 AU. With a 4 yr long mission, the 
three planets are identified, but their masses and SMAs 
are poorly constrained (relative errors are approximately 
10% for the SMA and 30% for the mass). There is some 
confusion between planets 1 and 2, and a correlated er- 
ror between the mass and SMA for each of the planets. 
These results show that identification of the potentially 
habitable inner planet in this case requires about 4 yr. 

Figure m shows how well the planets masses and SMAs 
are constrained for combined coronagraphy -I- astrometry 
missions with 1 yr, 2 yr and 4 yr durations. After only 1 
yr, the three planets are clearly identified, even if for the 
two inner planets, there is no strictly positive lower limit 
to the planet mass (detection is then confirmed by the 
coronagraphic images). The inner planet is unambigu- 
ously identified as a low-mass (< 5M Earth) planet in the 
star's habitable zone. The simulation therefore demon- 
strates that a 1 yr mission duration (with measurements 
at only six epochs) is sufficient to identify high-priority 
targets for follow-up. If only astrometric measurements 
were used, it would take approximately 4 yr to reach the 
same conclusion. This figure also shows that the outer 
planet's mass is relatively well constrained even though 
it has only covered about one quarter of its orbit. With 
a 2 yr mission duration, the relative errors in the masses 
and SMA estimates are better than what is obtained in 
the astrometry only mission in 4 yr. 

3.4. Discussion 

Results shown in Figure [6] demonstrate that a com- 
bined astrometry -I- coronagraphy mission can be de- 
ployed with no a priori knowledge of the target locations, 
and quickly (in 1-2 yr) identify potentially habitable 
planets. This is a significant advantage over the astrom- 
etry only case for which, as described in the Appendix, 
proper identification of planets requires that all planets of 



the system are either identified by RV or cover > 80% of 
their orbit during the measurement span. This fast iden- 
tification would allow the first phase of the an imaging 
mission -f- astrometry mission to perform high-cadence 
observations of a large number of targets, while the sec- 
ond phase of the mission could focus on acquiring high- 
quality spectra and astrometry of the low-mass planets 
identified in the habitable zone of the sample stars. 

Figure |6] illustrates that a key advantage of the com- 
bined measurement, beyond the improvement in mass 
and orbital parameters estimations, is the very low prob- 
ability of obtaining a false solution when solving for the 
orbital elements and masses of multiple planets. With 
astrometry alone (Figure [5]) or combined RV and as- 
trometry measurements (see the Appendix), the prob- 
ability of obtaining false solutions is non-negligible and 
rapidly increases as the measurement S/N or number of 
measurements decrease. The combined astrometry and 
coronagraphy measurement is very robust against such 
errors even for low astrometry S/N. This good perfor- 
mance in determining the masses and orbits of planets in 
multiple systems is particularly relevant i n light o f recent 
KeplerSpaceTelescope results. Lissau er et al.l (|2011l ) 
find that 408 of the 1200 recently released Kepler can- 
didate planets are in multiple transiting systems. Given 
that the orbits of planets in these systems must be co- 
aligned to « 1° to be detected by Kepler, this suggests 
that all or nearly all planets may be in multiple systems 
if the orbital inclinations of planets are aligned to only 
w 10° in each system. 

While the planets chosen for this example are expected 
to have similar reflected light contrasts (the increased 
planet mass as a function of angular separation results in 
roughly constant surface area divided by squared SMA) , 
the astrometric signal is dominated by the outer planet 
(32 times stronger signal than the inner planet). The re- 
sults presented in this section thus indicate that the rel- 
atively small astrometric signature of a habitable planet 
can be reliably separated from the much stronger signal 
of outer planets. The ability to identify and measure the 
mass of planets as a function of planet mass and separa- 
tion is expected to follow the trends discussed in Section 
12.61 at the small separations, the coronagraph's IWA lim- 
its sensitivity, while the coronagraph contrast, and possi- 
bly mission duration, limit sensitivity for outer planets. 
With multiple planets, confusion between planets may 
also become an issue, and could be especially serious for 
a system viewed edge-on, for which the number of in- 
dependent measurements is essentially halved. Planet 
colors and brightness as a function of illumination phase 
(ignored in this paper) would then need to be taken into 
account to help identify individual planets. We note that 
the proposed astrometry -I- coronagraphy measurement 
scheme is however well suited to this challenging configu- 
ration, compared to a much more degenerate astrometry 
only measurement scheme. 

We note that our simulations do not assume avail- 
ability of RV measurements, which have been shown to 
greatly mitigate confusion issues for an astrometry-only 
mission (see the Appendix). RV measurements are also 
expected to be valuable for an astrometry -|-coronagraphy 
mission, as massive distant planets would not be visible 
in reflected light in the coronagraphic images, but would 
produce large- amplitude long-period astrometric signals. 
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Unknown planets within the coronagraph IWA will pro- 
duce astrometric noise for the combined astrometry + 
coronagraphy approach described in this paper. These 
planets will have shorter period, and, at constant mass, 
a smaller astrometric signature than planets outside the 
IWA. They can therefore only compromise the solution 
if significantly more massive than the planets for which 
the mass is to be measured. This problem can be mit- 
igated by a larger number of measurements, especially 
at high cadence, to (1) identify the inner planet (s) from 
astrometry alone and (2) average down the astrometric 
noise from unknown short-period planet(s). RV mea- 
surements, if available, would be very well suited to help 
identify such planets. 

4. DEEP WIDE-FIELD IMAGING SENSITIVITY WITH A 
DIFFRACTIVE PUPIL TELESCOPE 

The effect of the dots on the telescope primary mirror 
on general-purpose wide- field imaging is eva luated in this 
section. The baseline design described in iGuvon et al.l 
([2012) is adopted in this paper: a 1.4 m telescope with 
dots covering 1% of the primary mirror area, with a 0.29 
deg^ field of view camera with 44 mas pixels observing in 
visible light. We assume in this section that this system 
is observing a Sun-like star at 6pc {my = 3.7). 

4.1. Light Lost to the Dots 

The light lost due to the dots is the sum of the light 
directly absorbed by the dots and the light they diffract 
out of the point spread function core into a wide halo 
of spikes. This second component is diffracted out to 
large distances, and is considered lost for science. Both 
quantities are equal to the fractional area of the pupil 
occupied by the dots. In the baseline design chosen, the 
total photometric loss is therefore 1%-|-1% = 2%. We 
note that this loss is small compared to the sensitivity 
gain offered by adopting an unobstructed pupil, as the 
photometric sensitivity loss due to spiders and central 
obstruction in an on-axis telescope design would be larger 
than 2%. 

4.2. Additional Background Due to Light of the Central 
Star Diffracted by the Dots 

The central star's diffraction spikes extend over most 
of the wide-field image. Although their contrast relative 
to the central star is faint (approximately 10~® along 
the spikes), the central star is much brighter than other 
sources in the field. 

Over 50% of the field, the additional diffracted light is 
less than 2.8 photons per day per pixel. Table |4] shows 
how, for a my = 3.7 central star, the light from the cen- 
tral star compares with the zodiacal light background. 
The 50 percentile line shows that the median diffracted 
light surface brightness is 0.03% of the zodiacal light 
brightness. Over 95% of the field, the additional light in- 
troduced by the dots on the primary mirror is less than 
1% of the zodiacal background. If the telescope were 
pointed on Sirius, 95% of the field sensitivity would be 
photon-limited from zodiacal light, while over 5% of the 
field, the photon noise from the spikes would limit sen- 
sitivity. The zodiacal light photon-noise-limited point 
source detection limit for the 1.4 m diameter telescope 
is mv ~ 32 for a 2 day unfiltered exposure, and would 



TABLE 4 
Scattered Light Due to Primary Mirror Dots 

Field Percentile^ Flux{photon.s~^ .pixel~^)^ Flux/Flux^o^ 



1 % 


8.28 lO"'^ 


7.66 10-5 


5% 


1.08 IQ-s 


9.96 10-5 


10% 


1.28 10-5 


1.19 10-* 


20 % 


1.68 10-5 


1.55 10-* 


50 % 


3.23 10-5 


2.99 10-* 


80 % 


9,27 10-5 


8.56 10-* 


90 % 


3.06 10-* 


2.82 10-3 


95 % 


1.19 10-3 


1.10 10-2 


99 % 


4.35 10-2 


4.02 10-1 



^ All fluxes are measured in a 3x3 arcmin box centered 6 arcinin 

from the optical axis. Scattered light becomes smaller at larger 

separations due to the Airy pattern contribution and the diffraction 

envelope of the spikes (equal to the diffraction pattern of a single 

dot). 

^ All values in this table assume an unfiltered exposure (zero point 

= 5.5^" photon. sec~^), and a my = 3.7 central star 

'^ Assuming mv = 22.5 magn.arcsec~ zodiacal light background 

level 

be mv ~ 34 for a full set of 32 observations, each 2 day 
long. 

The diffraction spikes can also be numerically removed 
from the wide-field image to a high accuracy, as they 
are, to first order, static on the detector. Since the ob- 
servation mode has the telescope slowly rolling during 
the observation, background sources are rotating on the 
detector while the spikes are static. A median of the 
images acquired during the roll contains only the spikes 
(and the zodiacal light background), and can then be 
subtracted to each individual frame prior to de-rotation 
and co-addition. 

5. CONCLUSIONS 

The DPT design allows simultaneous deep wide-field 
imaging, coronagraphic imaging, and astrometric mass 
determination of exoplanets around nearby stars. We 
have shown and quantified in both this paper and our 
previous publication that the measurements can be per- 
formed without impacting each other's performance: the 
coronagraphic and deep wide-field observations are not 
significantly affected by the presence of dots on the pri- 
mary mirror. The potential cost and complexity saving 
offered by combining three separate missions into one 
observatory could be significant, and we are thus cur- 
rently evaluating the feasibility and performance of this 
concept through both simulations and laboratory demon- 
strations. 

We note that the benefits of performing simultane- 
ously the three observations extend beyond the quantifi- 
able metrics (planet mass and orbital parameters uncer- 
tainties) discussed in this paper. Together, astrometric 
measurements and coronagraphic images are a powerful 
way to avoid confusion problems in complex planetary 
systems, where a large number of planets, comets, dust 
clouds may be present. Astrometric measurements may 
also allow recovery of planets which are at the limit of de- 
tection in the coronagraphic images, for example because 
of residual structure in the scattered light halo: the astro- 
metric measurements could confirm or eliminate features 
which could be either a planet or a speckle. This is par- 
ticularly beneficial if the astrometric S/N is high but the 
contrast of the planet makes its identification in corona- 
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graphic images ambiguous - images acquired at different 
epochs can then be co-added following the planet loca- 
tion to reach deeper contrast (by the square root of the 
number of measurements) than otherwise possible. This 
scheme becomes especially valuable for planets at the 
outer edge of the habitable zone, which are fainter but 
have stronger astromctric signatures than planets closer 
in. 

We note that these benefits do not require strict si- 
multaneity of the astrometric and coronagraphic mea- 
surements as long as the two types of measurements can 
be connected together. Separate astromctry and direct 
imaging missions, if operating within a few years, would 
thus offer the same advantages, provided that the as- 
tromctry measurements are not performed after the di- 
rect imaging measurements (this would not allow prior- 
itization of targets for follow-up spectroscopy of poten- 
tially habitable Earth-mass planets with the direct imag- 
ing mission). 

Acquiring deep wide-field images around nearby stars 
may also reveal extended debris disks (analog to the 
Kuiper Belt, scattered disk, or Oort Cloud in the so- 
lar system), and therefore advancing our understanding 
of the planetary system architecture and history. Fi- 
nally, we note that the spikes image in the wide-field 
camera are a low-resolution spectrum of the central star. 



which, when acquired simultaneously with the corona- 
graphic and astrometric measurements, may allow addi- 
tional calibration. 

Our study also shows that coronagraphic imaging and 
astrometry allow stellar mass measurements to percent- 
level precision for any star around which planets can 
be identified, allowing calibration of the stellar mass- 
luminosity relationship. 

A laboratory effort is under way to validate the concept 
and demonstrate that sub-/ias level astrometric precision 
is attainable with a wide-field imaging telescope. A sep- 
arate effort is aimed at laboratory demonstration that 
the DTP concept is compatible with high contrast coro- 
nagraphic imaging. Results from these experiments will 
be reported in future publications. 
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APPENDIX 
FINDING PLANETS FROM ASTROMETRY ONLY 

In this Appendix, we briefly review what can be detected by astrometry alone, a scenario that has been studied as 
part of the Space Interferometry Mission (SIM) project. In addition to astrometry, it is here assumed that all the 
nearby stars will have been monitored by RV measurements. While RV cannot detect 1 Mearth mass planets in a 
1 yr orbit, due to stellar RV noise, RV because of its long time baseline will have found larger planets Jupiters and 
Neptunes with long orbital periods. A detailed study of what is possible when 5 yr of precision astrometric data plus 
wl5 yr of 1 m.s^^ RV data was conducted during a "double blind" simulation of a 5 yr SIM. We briefly summarize in 
this section t he lessons lea r ned d u ring th i s doub le blind study. A more detailed description of the double blind study 
is provided in lTraub et all pOlol ): iTraubl pOlol ). 

The first conclusion is that astrometric detection of exo-Earths requires a mission S/N=6. Here S/N is defined 
as signal-to-noise where signal is the semi-major axis of the astrometric orbit and noise is the Ict rms error of the 
whole 5 yr data set. If there were 100 measurements during the mission and 12 parameters are solved for (five stellar 
parameters: position, proper motion, and parallax; six orbital parameters and the mass of the planet), the Ict 5 yr 
mission noise is the single epoch noise divided by VlOO — 12. With an S/N> 6, the false alarm probability (FAP) is 
< 1%. FAP of 1% means that if there were no planets, and 100 stars are searched, noise will mistakenly be identified 
as a planet once. Second, at S/N= 6, the mass of a 1 Mearth planet would have a ~25% Ict error. The orbital phase 
in radians has an error that is «0.5 of the mass. The period is measured more accurately for short-period planets, a 
planet with a 2x shorter period have its period measured twice as accurately. Last of all, the ability to find a planet 
and measure its mass was independent of the planet's period as long as the period was less than «80% of the mission 
length. The sensitivity degraded with longer periods, by about 2x when the planet period equaled the mission length. 

The main purpose of the double blind test was to examine the astrometric detection of multiple planets. In general 
the presence of multiple planets does not present a problem for astrometric detection if the planet orbital frequencies 
are separated more than 1/mission length. The double blind test used a large variety of multiple planet systems and 
in general confusion in the frequency domain was not an issue except for planets whose periods were longer than the 
mission length. The parallax effect has a 1 yr period and planets whose orbital periods were between 0.9 and l.lyr 
would have part of their motion absorbed into the parallax solution. Again out of a randomly selected 100 planetary 
systems only one had an earth like planet close enought to l.Oyr to present a problem but by chance its orbital plane 
was almost orthogonal to the parallax signature and its mass was properly measured. If a multiple planet system had 
one large outer planet, such as Jupiter in a 12 yr orbit, it was still possible to "fit" for the 12 yr planet and after that 
find all the planets with shorter periods whose S/N was > 6. The problem arose when there were two outer planets, 
for example with 12 and 20 yr periods. In this situation fitting a single long-period planet to the data when there were 
actually two long-period planets could result in significant residuals. The fitting error in general is an "arc" of motion 
over the 5 yr data period. But because proper motion is also fitted at the same time, the arc becomes a closed curve, 
mimicking a planet with a period slightly less than the mission lifetime. With experience the teams participating in the 
double blind study were able to recognize this about 70% of the time and were then able to recover both long-period 
planets and the shorter period planets. When discussing the ability of astrometry to detect a planet in a multiplanet 
system, there are two statistical metrics, one is called completeness, what fraction of the planets whose S/N> 6 and 
had periods less than 5 yr were actually found and the second is confidence, what fraction of "claimed" discoveries 
were actual planets. The result of the double blind test was very encouraging in that both completeness and confidence 
were above 90%. 



